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Abstract 

The fundamental BCJ-relation is a linear relation between primitive tree amplitudes with 
different cyclic orderings. The cyclic orderings differ by the insertion place of one gluon. 
The coefficients of the fundamental BCJ-relation are linear in the Lorentz invariants 2 pipj. 
The BCJ-relations are well established for pure gluonic amplitudes as well as for amplitudes 
in = 4 super-Yang-Mills theory. Recently, it has been conjectured that the BCJ-relations 
hold also for QCD amplitudes. In this paper we give a proof of this conjecture. The proof is 
valid for massless and massive quarks. 


1 Introduction 


Amplitudes in QCD are often computed by a decomposition into a sum of smaller pieces, called 
primitive amplitudes [Hl|2]|. The primitive amplitudes are gauge invariant, colour-stripped and 
have a fixed ordering of the external legs. Non-trivial relations between different primitive tree 
amplitudes are a fascinating topic and have important applications. For pure gluonic primitive 
tree amplitudes these relations are by now well-studied. Invariance under cyclic permutations 
is trivial. The first non-trivial relations are the Kleiss-Kuijf relations 0, which follow from the 
anti-symmetry of the colour-stripped vertices. More interesting are the Bern-Carrasco-Johansson 
relations (BCJ-relations) 0. The full set of the BCJ-relations follows from the so-called funda¬ 
mental BCJ-relations Q. The fundamental BCJ-relation is a linear relation between primitive 
tree amplitudes with different cyclic orderings. The cyclic orderings differ by the insertion place 
of one gluon. In the fundamental BCJ-relation the coefficients of the relation are linear in the 
Lorentz invariants 2 PiPj. The BCJ-relations are known to hold for pure gluonic tree amplitudes 
and for tree amplitudes in = 4 super-Yang-Mills theory. The BCJ relations have been proven 
first with methods from string theory [[610 and later within quantum field theory with the help 
of on-shell recursion relations ©El®. On-shell recursion relations require a certain fall-off be¬ 
haviour for large momentum deformations. For amplitudes in = 4 SYM it is essential that 
not only the (bosonic) momentum components but also the (fermionic) Grassmann components 
are shifted. The required fall-off behaviour has been established for pure gluonic tree amplitudes 
and amplitudes in = 4 SYM in [fT0l - IT2l . Furthermore BCJ relations have been derived for a 
pair of massive scalars and an arbitrary number of gluons [|T3l . 

It is natural to consider primitive tree amplitudes in QCD as well, involving in addition to 
gluons massless and/or massive quarks ([141115 1. The fundamental BCJ-relation singles out three 
particles, which we will label 1, 2 and n. In the fundamental BCJ-relation the positions of 
particles 1 and n are fixed, as there the positions of the remaining particles 3 to (n — 1). In the 
cyclic order particles 1 and n are adjacent. In the fundamental BCJ-relation particle 2 is inserted 
in all possible places in the cyclic order between 1 and n, but not between n and 1. Recently, 
Johansson and Ochirov conjectured lfl4l that the fundamental BCJ-relations hold for primitive 
tree amplitudes in full QCD as well, provided particle 2 is a gluon. In this paper we prove this 
conjecture. The major part of our arguments is not specific to four space-time dimensions. Only 
in the explicit definitions of the momentum shifts we use spinor representations corresponding to 
four space-time dimensions. With a suitable generalisation of these momentum shifts our proof 
will carry over to D space-time dimensions. 

This paper is organised as follows: In section © we introduce the conjecture, which we 
are going to prove and give an outline of the proof. Section © carefully defines momentum 
deformations through three-particle shifts. This is a necessary technical preparation for our proof. 
In section © we discuss the large z-behaviour of the deformed fundamental BCJ-relation under 
the three-particle shifts and show that there is no contribution from infinity in BCFW-recursion 
relations. In section © we prove the fundamental BCJ-relation by induction. Our conclusions 
are given in section ©. In an appendix we collected some technical details on certain three- 
particle shifts with massive quarks (appendix [A]) and a proof on the large z-behaviour in the 
eikonal approximation (appendix [B]). For the convenience of the reader we also included the 
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cyclic-ordered Feynman rules relevant to primitive amplitudes in appendix O 


2 Overview 

Tree amplitudes in QCD with any number of external quarks can be decomposed systematically 
into gauge-invariant primitive amplitudes with a fixed cyclic ordering []2][m|T7]|. Let us consider 
a tree-level primitive QCD amplitude with n external particles, out of which n q particles are 
quarks, n q particles are anti-quarks and n g particles are gluons. We have the obvious relation 

n — n g + 2n q . (1) 

Without loss of generality we may assume that all quarks have different flavours. The quarks may 
be massless or massive. In this paper we consider amplitudes with at least one gluon (n g > 0). 
This excludes the case, where all external particles are either quarks or anti-quarks (n g — 0 and 
hence n = 2n q ). We are mainly interested in the case, where there is in addition to gluons at least 
one quark-anti-quark pair ( n q > 0). The pure gluonic case (n q — 0) is not excluded, but already 
well studied. The tree-level primitive QCD amplitudes are cyclic-ordered. We denote such an 
amplitude with the cyclic-order (1,2,by 

A n (1,2,...,«). (2) 

The amplitudes satisfy the Kleiss-Kuijf relations [|3|. In order to state the Kleiss-Kuijf relation 
let 

a = (ai, ...,cc/), p = (pr, ...,(3 n _ 2 -y) (3) 

and fr = (P„_ 2 ..., Pi). The Kleiss-Kuijf relations read 

A, ; (l,p,2,dc) = (- l) n ~ 2 ~ j £ A„( 1,2,0!,..., 0 ,^ 2 ). (4) 

cealll j3 r 

Here, a III p r denotes the set of all shuffles of a with p 7 , i.e. the set of all permutations of the 
elements of a and p 7 , which preserve the relative order of the elements of a and of the elements 
of p r . A special case of the Kleiss-Kuijf relations is the situation, where the set P contains only 
one element. In this case the Kleiss-Kuijf relation reduces to the U (1)-decoupling identity 

(CTl, 02: •••; 1 ? ^) 0; (5) 

i 

where the sum is over the cyclic permutations of the first (n — 1) arguments. The Kleiss-Kuijf 
relations in eq. © allow us to fix two legs at specified positions. We will use this freedom to fix 
a particle at position 1 and a second particle at position n. Since we also assume that there is at 
least one gluon, let us label this gluon by 2 g . In | fl4| Johansson and Ochirov conjecture that 

n— 1 / n \ 

£ £ 2p2Pj]A n (l,3,...,i,2 g ,i+l,...,n-l,ri) = 0. (6) 

i =2 \j=i +1 / 
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In this paper we prove this conjecture. An equivalent formulation of eq. © is 

n— 1 

i=2 

which follows from eq. © by momentum conservation. Eq. © is the fundamental BCJ-relation 
013 for tree-level primitive QCD amplitudes. It is well known that eq. © holds in the pure 
gluonic case. It is also know that eq. © holds for amplitudes with one quark-anti-quark pair 
(n q — 1) in massless QCD. This follows from the fact that these amplitudes are identical to the 
corresponding amplitudes in = 4 super-Yang-Mills theory, and the fact that the BCJ-relations 
hold for the latter liSl ffT 1 0. The aim of this paper is to show that eq. © holds more generally in 
(massless or massive) QCD. 

In order to prove the fundamental BCJ-relation we will make use of on-shell recursion re¬ 
lations (or Britto-Cachazo-Feng-Witten recursion relations, BCFW-recursion relations for short) 
[f2T1 . Within on-shell recursion relations some distinguished external momenta are deformed, 
such that momentum conservation and the on-shell conditions are respected. Fet us consider a 
deformation of the momenta pi, P 2 and p n depending on a variable z. This is called a three- 
particle BCFW-shift m. Since for our problem there are three distinguished particles 1, 2 g and 
n, a three-particle shift is more natural as compared to a more conventional two-particle shift. It 
will turn out that a three-particle shift will simplify the proof. We denote the deformed momenta 
by 



Y, 2 P 2 Pj ) A n (1,3, 2 g ,i+ l,ri) = 0, 


( 7 ) 


V7=l 


Pl(z), p 2 (z), p n {z). 


(S) 


We further require 

M°)=Ph P2{0)=P2 , Pn(0)=Pn■ (9) 

For j 4 1.2.// we simply set P jiz) — Pj. We introduce the quantity 

I«{z) = if I 2 P 2 PJ 
i=2 \j=i +1 

For z — 0 the expression I„(z) reduces to the left-hand-side of eq. ©. I n (z) is clearly a rational 
function of z. We have to show 


I A n (l,3,...,i,2 g ,i+l,...,n-l,h ). 


( 10 ) 


4(0) = 0, 


( 11 ) 


or equivalently 



z=0 


( 12 ) 


We do this in two steps: We first show that I n (z) has no pole at z — °°, or equivalently, that I n {z) 
falls off for large z at least with 1/z. In the second step we use induction and BCFW-recursion in 
order to prove eq. (fTTI) . 

'it is worth noting that all tree amplitudes in massless QCD can be obtained from tree amplitudes in ZK[ = A 
super-Yang-Mills theory Ifl8ll20l . 
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3 Momentum deformation 


In this section we review the three-particle BCFW-shift for the external particles 1, 2 g and n. 
The types of particles 1 and n may either be a quark, an anti-quark or a gluon. However, particle 
2 g is required to be a gluon. Quarks and anti-quarks may be massive or massless. Therefore 
we will treat the general case with arbitrary masses for particles 1 and n. The case of massless 
quarks is included as the special case m\ — m n — 0. This section may appear at a first reading a 
little bit technical, but actually it will pay off: We define the momentum deformations in such a 
way, that we obtain the same large z-behaviour independently of the helicity configuration and 
independently of the masses of the particles. Most of the technical details are related to massive 
quarks and it might be advantageous at a first reading to pay attention to the massless case only. 
In the massless case the formulae simplify considerably. 


3.1 Spinor definitions 

For the definition of massive spinors we follow lf23l . Let us consider two independent Weyl 
spinors |g+) and (q + |. These two Weyl spinors define a light-like four-vector 


f\q+)- (13) 

This four-vector can be used to associate to any not necessarily light-like four-vector p a light¬ 
like four-vector p K . 


P 


b 



(14) 


The four-vector p v satisfies (p b ) 2 = 0. Note that p° depends implicitly on \q+) and (q+ |. The 
two Weyl spinors \q+) and (q-\- | are used as reference spinors in the definition of the polar¬ 
isations of the external particles. For quarks with momentum p we take the spinors u and u 
as 


u(-) = 

M (+) = 

The spinors v and v 

v(-) = 
v(+) = 


(p b + \q ) ^ 

(p b - \q+) ^ 

are given by 


+ m) \q~), 
+ m) |<?+), 


1 

(p b + \q~) 

1 

(p b - k+) 


Q/-m) | q+), 


«(“) = , , I b , (<? + I (// + m ) • 

(■ q+\p-) 


V(+) 

v(-) 


(q-\p^+) ^ 

1 

{q + lp 9 -) ^ + 


I (7/ - m), 
I (7 1 ~ m) . 


(15) 


(16) 


We label the helicities as if all particles were outgoing. As a consequence, the spinors u(k) 
and v(X), which correspond to particles with incoming momentum, have a reversed helicity 
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assignment. In the massless limit the definition reduces to 


u(~) = v(-) = \p+), u{+)=v{+) = {p+ |, 

«(+) = v(+) = \p~), *(-) =v(-) = {p-\, (17) 


and the spinors are independent of the reference spinors | q+) and (q+ |. 

For massive fermions the reference spinors are related to the quantisation axis of the spin 
for this fermion, and the individual amplitudes with label + or — will therefore depend on the 
reference spinors |<?+) and {q+ |. It should be noted that the spinors for the massive fermions 
depend both on |g+) and (q + \ : For the spinors with helicity “+” there is an explicit dependence 
on | q+), while the dependence on (q+ | enters implicitly through p’. In a similar way we find 
that the spinors with helicity ” have an explicit dependence on (q+ |, while the dependence 
on |</+) enters implicitly through p ’. 

It is easy to relate helicity amplitudes of massive quarks corresponding to one choice of 
reference spinors to another set of reference spinors. If \q+) and (q + is a second pair of 
reference spinors we have the following transformation law 


f u(+,q) \ = ( cn c 12 \ f u(+ , q) \ 

V / V C 21 c 22 ) V «(->«) ) ’ 


where 


Cll 


(g- W\q~) ^ = m (m) 

(qp b )\p b q] ’ CU (#)(A)’ 


m\qq\ ^ = (q + \l/\q+) 
qp h ] [p J q]' C22 [qp^pS) ' 


( 18 ) 


(19) 


Here, p’ denotes the projection onto a light-like four-vector with respect to the reference vector 
\(q + | , y t '|<3 : +). Similar, we have for an amplitude with an incoming massive quark 


f «(+,<?) \ = / cn -cn \ ( u(+,q) \ 

V ) v _C 21 c 22 ) V u (~iq) )' 


( 20 ) 


Similar formulae exist for the spinors v and v and can be obtained by the substitution u — > v, 
u — y v and m —> —m. 

For the polarisation vectors of a gluon with momentum p we take 

£+ = (p+\Uq+) £ = (q+\j M \P+) (21) 

11 V2(q-\p+)' M y/2{p + \q~) 

The dependence on the reference spinors |< 7 +) and (q+ |, which enters through the gluon polar¬ 
isation vectors will drop out in gauge invariant quantities. 


3.2 Decomposition of massive four-vectors into light-like four-vectors 

The external momenta of particles 1 and n may be massive or massless. In the case where they 
are massive (either one of them or both) we would like to write them as a linear combination of 
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two light-like four-vectors l\ and The two light-like four-vectors l\ and l n are constructed as 
follows [1241 1251: If p\ and p n are massless, l\ and l„ are given by 


— Pi j hi — Pn- 


If pi is massless, but p n is massive one has 


h=Ph hi = -tt\P\ +Pn, OCl 



2p\Pn 


The inverse formula is given by 


( 22 ) 


(23) 


Pi=h, p n = aih+hi- 
If both p\ and p n are massive, one has 

h = - - \~^~{pi -ttnPn) , hi = 7 --- (-UlPl+Pn)- 


1 (X. | (X,; 


1 — CXiCt,, 


(24) 


(25) 


aj and a„ are given by 


2pip„-sign(2pi;^)v / A _ 2pip n - sign(2p 1 p„)v / A 
CXl 1 5 ct n 


2 P\ 


2 p 


(26) 


Here, 


A = {2pip n f - Ap\p 2 n . (27) 

The signs are chosen in such away that the massless limit p\ —y 0 (or p\ —> 0) is approached 
smoothly. The inverse formulae are given by 

Pi — ll T Ct/ihii Pn CCt/i T In- (28) 

The light-like four-vectors l\ and /„ define massless spinors |/i+), (/| + |, |/„+) and (/„+ |. 

3.3 On the choice of the reference spinors 

Particles 1 and n may be massive quarks or anti-quarks and we have to make a choice for the 
reference spinors. In the massless case, the amplitude will be independent of the choice of the 
reference spinors and this section is of no further relevance. However, if particle 1 (or particle 
n) is massive, the choice of the reference spinors will define the spin quantisation axis and the 
amplitude will depend on this choice. It is always possible to convert to a different basis with the 
help of eqs. (fT8T)-(l20l). 

In section (13.21) we have constructed the spinors |/i+), | l n +), {h + \ and (l n + |. For generic 
momenta p\ and p n . the ket-spinors |/i+) and |/„+) span the two-dimensional space of holo- 
morphic spinors (or ket-plus-spinors). Similarly, the bra-spinors (/i + | and (/„ + | span the space 
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of anti-holomorphic spinors (or bra-plus-spinors). For particle 1 we parametrise the reference 
spinors |<?i+) and {q\ + | as 

|#t+) — K«+) + ^-t |h+) > (#i + | — (^n+| + ^t (h+| j (29) 

with one parameter Xj. For particle n we parametrise the reference spinors | q„+) and (q„ + | as 


|#77+) — 1^1+) + hn |/«+) , (#71+1 — (h + | +^77 (^7I+| ) 

with one further parameter X n . Working out | p\+), \p\+), (p\ + | and (p^ + | one finds 


(30) 


p\+ 


- c\ (|/t+) — oc„A.i |/,7+)), 

Pn + ^ = C tl (1^77 + ) — OCl A.,7 |/l+)) , 


with 


c l = 


1 + 0C ;I A, ^ 


p\ + 

pi+ 


Cn — 


— Cl ((/l + | — (^71+|) , 

= C^((/fj- 1 -| 0^1^77 (^1 + 1 ) 5 

1 




(31) 


(32) 


The variables ai and a n have been defined in eq. (l26l) . We will use the freedom to choose A.| and 
X n to compensate a restriction on the possible BCFW-shifts for massive particles. The shifted 
spinors of the massive particles have to satisfy the Dirac equation, as well as orthogonality and 
completeness relations. This restricts the z-dependent part to be proportional to the correspond¬ 
ing reference spinors f[23l . This means if we shift a massive spinor 

«(-) = w(-)-z|'n+), (33) 

the spinor |r|-|-) has to be proportional to |#+): 

|r|+) = k |#+). (34) 

Similarly, if we shift 

m(+) = m(+) + z (t|+| , (35) 

we have to require that 

(tl+| = K (#+|- (36) 


3.4 BCFW-shifts 

We now provide the explicit formulae for the three-particle shifts. We have to consider all pos¬ 
sible helicity configurations for the three particles 1, 2 g and n. In all cases the deformations are 
defined in such a way, that the external polarisations spinors and vectors give the best possible 
large /-behaviour. This is a z -3 -behaviour if the three particles 1, 2 g and n are all gluons, a 
z -2 -behaviour if one of them is a quark or an anti-quark and a z -1 -behaviour if two of them are 
quarks or anti-quarks. As particle 2 g is required to be a gluon, the case where all three particles 
are quarks or anti-quarks is not possible. The large z-behaviour of the external polarisations is 
summarised in table CD). 
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particles 1, 2 g and n 

large z-behaviour 

3 gluons 

z” j 

2 gluons, 1 quark/anti-quark 

z " 2 

1 gluon, 2 quarks/anti-quarks 

z - 1 


Table 1: The large z-behaviour of the external polarisations under the three-particle shifts. 


3.4.1 The helicity configuration 1 + , 2+, n 

For the helicity configuration l + ,2we shift u i(—), « 2 (“) and while mi(+), « 2 (+) 

and u n {—) remain unchanged: 

ui(-) = ui(^)-zyi\p b n +), Un(+) = Un(+)+zyi{p\ + \+zy2{P2 + \, 

u 2 {~) = u 2 (-)-zy 2 \p b n +), (37) 


where y | and y 2 are two non-zero constants. For massive particles we have to require, that the 
shift is proportional to the corresponding reference spinors. Therefore we have to require that 
the system of equations 


Ki |^l+) 


yi 



K/i (< 7 » T | — y i 


p\ + 


+y 2 (P2+I, 


(38) 


has a solution for some non-zero constants Ki and K 2 . In appendix (IA.1I) we show that this is 

the case. The spinors u \{—) and wi(+) correspond to an on-shell particle with mass m\ and 

four-momentum 

f\ = p^-\zy\(p\ + \f\pi+) ■ (39) 

The spinors u 2 {—) and u 2 (+) correspond to an on-shell particle with zero mass and four-momentum 

P 2 = P 2 ~\zy2( K P 2 + \f\Pn+)- (40) 

The spinors u n (—) and u n {+) correspond to an on-shell particle with mass m„ and four-momentum 

K = P% + ^m(p\ + \ r f\p l ’n+) + ^zy 2 (p 2 + \f\p b n +y (41) 


3.4.2 The helicity configuration 1 + ,2 g ,n 

For the helicity configuration l + ,2 g ,n~ we shift u i(—), u 2 (+) and u n (+), while u i(+), u 2 {—) 
and u n (—) remain unchanged: 

«i(-) = «t(-) ~zy 2 \p 2 +) —zy n \Pn+), « 2 (+) = « 2 (+) +zy 2 (p\ + I, 

Mrc(+) = Mm(+) T Zyn{p\ + |; 
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(42) 









where y 2 and y n are two non-zero constants. For massive particles we have to require that the 
system of equations 


Kl |#1+) 


yi \p2+)+y n 



K/J (^7/7 ~F 



(43) 


has a solution for some non-zero constants Ki and K„. In appendix (IA.2I) we show that this is 
the case. The spinors u \{—) and ni(+) correspond to an on-shell particle with mass m\ and 
four-momentum 


*/J 

Pi 


Py-^zyi (p\ + \f\P2+} 



(44) 


The spinors U 2 (—) and m (+) correspond to an on-shell particle with zero mass and four-momentum 


P 2 = rt + \zy2(p\ + \f\p2+). (45) 

The spinors u n (—) and w„(+) correspond to an on-shell particle with mass m n and four-momentum 

K = p’ J n + ^zyn^p b i + \'f\Pn+) ■ (46) 


3.4.3 The helicity configuration 1 + , 2 g ,n + 

For the helicity configuration l + ,2 g ,n + we shift u \(—), M 2 (+) and u ,,(—), while u\(+), U 2 {—) 
and u n (+) remain unchanged: 

Mi(-) = Mi(-) -zyi\P2+), M 2 (+) = U 2 (+) +zyi{p\ + \+zy n (p b n + l 

Mn(-) = U„(-)-Zyn\P2+), (47) 

where yi and y n are two non-zero constants. For massive particles we choose 

l#l+) = k/7 + ) = |/?2+), (#1+| = (#n + | = {P2 + \ (48) 

as reference spinors. The spinors u \{—) and wi(+) correspond to an on-shell particle with mass 
mi and four-momentum 

fi = P^-\zy,(p\ + \r\P 2 +). (49) 

The spinors m 2 (—) and t< 2 (+) correspond to an on-shell particle with zero mass and four-momentum 

P2 = P%+^w{Pl + \i t \P2+) + ^Zyn(pn + \'f\P2+y ( 50 ) 

The spinors u „(—) and u n {+) correspond to an on-shell particle with mass m n and four-momentum 

K = ri-\zy n (pi + \f\p2+y ( 51 ) 
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3.4.4 The helicity configuration 1 + , 2+, n + 

For the helicity configuration l + ,2+,n+ we shift u i(—), u 2 (—) and u n (—), while u i(+), u 2 (+) 
and w„(+) remain unchanged: 

Ml(-) = Ml(-)-Z P2 pS h + ), 

U 2 {~) = u 2 (-)-z P b n p\ |ri+), 

Un{~) = U„(-)-Z p\p2 h+). (52) 

Here, |r(+) is an arbitrary spinor. For massive particles we choose 

kl+> = k«+) — l r l + ); kl + l — (<ln +1 = (tl + | (53) 

as reference spinors. The spinors «i(—) and «i(+) correspond to an on-shell particle with mass 
mi and four-momentum 

Pi = P\~\ z P2Pn (pi + IY'I t 1+}- (54) 

The spinors U2{~) an d u 2 (+) correspond to an on-shell particle with zero mass and four-momentum 

P 2 = Pi~\ z PnP\ (P2 + \f\n+)- (55) 

The spinors u n {—) and u n (+) correspond to an on-shell particle with mass m n and four-momentum 

Pn = Pn~\ Z P b lP2 (p b n + \f\ r [ + )- (56) 

Momentum conservation is satisfied due to the Schouten identity. 

3.4.5 The remaining helicity configurations 

The shifts for the helicity configurations 

(1- 2" «+),(l- 2+ «+),(l- 2+ 0,(1“ 2J,0 (57) 

can be obtained from the helicity configurations 

(l+,2+ n-),(l + ,2-,n-),(l+,2-,n+),(l+,2+ n + ) (58) 

by exchanging holomorphic and anti-holomorphic spinors. 
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4 Large-z behaviour 

We consider f„(z) for n > 4. I n (z) is a rational function in z. We have to show that I n (z) falls off 
at z — °° at least with l/z. We will distinguish the cases, where the three particles 1, 2„ and n are 

(i) three gluons, 

(ii) two gluons and one quark/anti-quark, 

(iii) one gluon and two quarks/anti-quarks, not belonging to the same fermion line or 

(iv) one gluon and a quark-anti-quark-pair belonging to the same fermion line. 

Let us recall the definition of I n (z): 

In(z) = t(± 2/3 2P 

7=2 \j=i+l 

We note that the factors (2p2pj) are at the worst linear in z. A sufficient condition is to show that 
each amplitude A„(l, 3, ...,i,2 g , i + 1, ...,n — l,n) falls off at z — °° at least with l/z 2 . We will 
show that this holds for the cases (i)-(iii). 

However, a l/z 2 -fall-off behaviour of the amplitudes is not a necessary condition. In fact, 
in the case (iv) the amplitudes fall off only with l/z. In this case we show through a more 
sophisticated argument that the full sum I n {z.) falls off at z — 00 with l/z. 

4.1 Three gluons 

Let us start with the case (l g: 2 g: n g ). The external polarisation vectors contribute a factor z -3 . 
The most critical contribution from the vertices and propagators comes from diagrams, where 
there are only three-gluon vertices along the z-flow. For these diagrams there will be along the 
z-flow always one more three-gluon vertex as there are propagators, giving a net factor of z 1 . 
Therefore we obtain from these diagrams a total contribution of z -3 -z = z -2 . If internally a 
gluon propagator is replaced by a quark propagator, we have to change at least two three-gluon 
vertices into quark-gluon vertices. This improves the estimate by a factor l/z. Similarly, the 
replacement of one three-gluon vertex by a four-gluon vertex results in an improvement in the 
z-behaviour by a factor l/z. We therefore conclude, that the amplitude falls off at z = °° at least 
with l/z 2 . 

4.2 Two gluons and one quark/anti-quark 

The arguments for the cases (l q / q ,2 g ,n g ) and (l g , 2 g ,n q / q ) are very similar to the three gluon 
case. Although the external polarisations contribute now only a factor z~ 2 . the estimate from the 
vertices and the propagators is now z°. Again, the worst diagrams are the ones with a maximal 
number of three-gluon vertices along the z-flow. However, in the case at hand we must have at 
least one quark-gluon-vertex along the z-flow, improving the estimate by a factor l/z. Again we 
see that the amplitude falls off at z = 00 at least with l/z 2 . 


I A n (i,3, ...,i,2 g ,i+ 1, ...,n - l,h). 


(59) 
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4.3 One gluon and two quarks/anti-quarks, not belonging to the same 
fermion line 

Let us now discuss the case (l q / q ,2 g ,n q i / q r) with one gluon and two quarks/anti-quarks, where 
the two fermions do not belong to the same fermion line. This sub-case is straightforward: Al¬ 
though the external polarisations contribute now only a factor z _1 , the estimate from the vertices 
and the propagators is now z 1 . This is due to the fact that we have to change yet another three- 
gluon vertex into a quark-gluon vertex. Again one concludes that the amplitude falls off at z — 00 
at least with 1/z 2 . 


4.4 One gluon and a quark-anti-quark-pair belonging to the same fermion 
line 


The case (l q ,2 g , n q ) and (\ q . 2 g , n q ), where the two fermions belong to the same fermion line, are 
more complicated. Power-counting gives now a factor z _1 from the external polarisations and a 
factor z° from the vertices and propagators. The individual amplitudes fall off as 1/z for large 
z. In this case we show, that the sum I n (z) falls off as 1/z for large z. The worst diagrams are 
the ones, where the z-flow of gluon 2 g goes only through three-gluon vertices before it couples 
to the quark line. We have to show that in the sum the leading z-behaviour of these diagrams 
actually vanishes. For the leading z-behaviour we can use an argument of Arkani-Hamed and 
Kaplan iTTOTl : For large z we may view particles 1, 2 g and n as highly energetic particles moving 
in a soft background. All vertices along the z-flow reduce in this limit to eikonal factors, except 
the one where the three branches of the z-flow meet. In order to see this let us start from particle 
2 and consider the first vertex particle 2 meets. This three-gluon vertex couples particle 2, a 
current containing only soft particles 

// oft = i; oft (k+i,...,/) ( 60 ) 

and a current containing the other hard particles 1 and n 

jj ard = y^ d (i,3,...,M+l,...,n). (61) 


In the Feynman rule for the three-gluon vertex we only have to keep the z-dependent terms, 
yielding for the cyclic order 2 g , 7 soft , 7 hard 


i [- (e 2 • 7 soft ) (p hard ■ Z 1 ^) + (p 2 • £ 2 ) (y soft • ,Z ard ) - 2 (e 2 • Z 1 ^) [p 2 ■ 7 soft ) 


(62) 


The contraction of p 2 with £ 2 vanishes: /3 2 • £ 2 = 0. Furthermore, the current Z iard is conserved 
and we have Phmd ■ Z 1 ^ = 0. This leaves the eikonal contribution 


_ yhard,amputated j 


2 P2-J 


soft 


(. P2 + Psoft) 


(63) 


with 


Ziard,amputated = / (p 2 + p soft ) 2 


(64) 
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We may then repeat the argument with the next three-gluon vertex. A similar argument can be 
given for the z-flow along the quark line. Let us start at particle 1 and let us assume that this 
particle is a quark. We consider the first vertex particle 1 meets. This is a quark-gluon vertex, 
connecting particle 1, a gluon current containing only soft particles 

7 s° ft = ^(3, (65) 


and a hard spinorial current containing the other two hard particles 2 and n: 

l/hard = yhard (jfc + 1 n - 1, h) . 

Let us further define the hard amputated spinorial current as 

yhard . (Pl+Psoft) + ™ y hard, amputated 

(P l + ll soft ) 2 —>n 2 

Again we may neglect soft momenta in the numerator and we find 

_^ yr_ P 1 m _yhard,amputated jsoft 

(pl+psoft ) 1 -™ 2 M 

2 pi -7 soft 


^ yhard, amputated _|_ u l(Pl m)”fV 

(P \ + P soft) 2 — m 2 1 , ,j \2 


/hard,amputated 
(pl+p soft) 2 - m 2 


( 66 ) 


(67) 


7 , s . oft . 


( 68 ) 


In the first term on the right-hand side we recognise an eikonal factor, the second term vanishes 
due to the Dirac equation. As before, we may repeat the argument with the next quark-gluon 
vertex. 

The argument for the branch with the external anti-quark at position n is identical and not 
repeated here. The (1 q, 2 g . «^)-casc is very similar and not discussed in detail. 

The eikonal factors go to a constant for large z and we are left with a quark-gluon vertex 
contracted for the (l ? ,2g,n^)-case with u.\, £2 and v n . Let us denote this contribution by 

O 3 = iuijftVn^. (69) 


The quantity O 3 falls off like 1 /z for large z. It is important to note, that O 3 occurs in every 
amplitude contributing to I n (z ) in the (l ? ,2^,n^)-case. It may therefore be taken out of the sum, 
and we have to show that the remaining sum goes to a constant for large z. The remaining 
sum involves only the Lorentz invariants 2 p 2 pj and the eikonal factors. The proof is given in 
appendix ©. 


5 The proof by induction 

In this section we prove the fundamental BCJ-relation by induction. With the preparations of 
section © and section © we can do this independently of the helicity configurations and the 
masses. This is possible, since we have for I n (z) for all helicity configurations and all masses 
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a 1/z-behaviour for large z. However, we would like to point out one subtle point for massive 
quarks: We would like to show that the fundamental BCJ-relations holds for all helicities of the 
massive quark. The naive way to show this would be to fix a spin quantisation axis through a 
choice of reference spinors |g+) and (q+ | and to show the BCJ-relation for the helicities “+” 
and ” with respect to these reference spinors. This is not what we are doing. The attentive 
reader of sections (13.4.1I) - (I3.4.41) might have noticed, that the and ’’-helicities refer to 
different reference spinors. This is o.k., since amplitudes with different spin quantisation axes 
are related through eq. (fT8l) and eq. (l20l) . Therefore it is sufficient to know two independent 
amplitudes (say “+”-helicity with respect to q and ’’-helicity with respect to q) in order to 
know all amplitudes with spin quantisation axes q and q. This remark applies to each external 
particle individually and covers all possible cases for the external particles 1 and n, where we 
can have out of these two particles either zero, one or two massive particles. In the latter case the 
masses may be equal or not. 

5.1 Induction start: The case n = 3 

To start the proof by induction we consider the case n — 3. Throughout this paper we work with 
complex external momenta. The external momenta satisfy momentum conservation 

P\+P2+P3 = 0, (70) 

and the on-shell conditions 

P\=m , p 2 = 0, p 3 = m . (71) 

Particle 2 g will always be a gluon and is therefore massless. Particles 1 and 3 may be massless 
or massive. In the massive case, particles 1 and 3 are necessarily a quark-anti-quark pair of 
the same flavour. Therefore particles 1 and 3 will have the same mass m. For n — 3 external 
particles the momentum configurations satisfying eq. (ITOl) and eq. (ITTI) are in general complex. 
The fundamental BCJ-relation reduces to 

2p 2 p 3 A 3 (l,2 8 ,3) = 0. (72) 

For generic external momenta A3 (1, 2 gl 3) is finite and 

2P2P3 = {P2 + Pi) 2 ~m 2 = p\~m 2 = 0. (73) 

5.2 The induction step 

We now show that 


0(0) = 0 


holds for j = n, provided it holds for all j < n. We start from eq. (fl2l) 



z=0 


(74) 


( 75 ) 
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where the contour is a small counter-clockwise circle around z — 0. Deforming the contour to a 
large circle at infinity and the residues at the finite poles z a ^ 0 we obtain 

/„(0) = S-^resf^) , (76) 

where B denotes the contribution from the large circle at infinity. In section © we have shown 
that I„(z) falls off at least with l/z for z— y °° and therefore 

5 = 0. (77) 

It will be convenient to introduce the following notation for the various factorisation channels: 

A n (1,2, ...,k,P\ -P,k+l,...,n- 1 ,h) = 

Y, A k+i (1,2 ,...,k,p) ~^A n _ k+] (- P,k+ \,...,n-\,h) , (78) 

together with the convention that the hatted quantities are evaluated at z — Za- The sum is over 

the helicity of the intermediate particle. Let us look at the z-momentum flow for a three-particle 
BCFW-shift. For each diagram we may divide the z-dcpcndcnt propagators into three segments. 
Each segment starts at the common vertex, where the z-dcpcndcnt momentum flow meets and 
goes outwards towards the particles 1, 2 g and n. We may use these segments to divide the finite 
residues into three groups and we write 


4(0) — 5i+52-|-5„, 


(79) 


with 

( Y, 2 p 2 Pj) i^A n (l,3,...,k,P\-P,k+l,...,i,2 g ,i + l,...,n-l,n), 

i=2 \j=i+l ) k =3 

n— 1 / n \ 

R 2 = E E 2 P2Pj 

i=2 \j=i +1 ) 

i n— 1 

x£ E A n (k + l,...,i,2 g ,i+l,...,l,P\-P,l+l,...,n-l,h,l,3,...,k), 

k= 2 l=i 

n —1 / n \ n—2 

R,i = E E 2 P2PJ) E A n(i,3,...,i,2 g ,i+l,...,k,P\-P,k+l,...,n-l,n). 

i=2 \j=i +1 / k=i 

(80) 


Let us first look at 5i. We may exchange the summation over i and k as 


n— 1 i 


E E /(*>*) 

i=2 k= 3 


n— 1 n— 1 


E E/ft*) ■ 

k =3 i=k 


(81) 
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One obtains 


Ri = ££( £ 2p 2 Pj)A n (l,3,...,k,P\-P,k+l,...,i,l g ,i+l,...,n-l,n ). (82) 

t=3 i=& V=/+l / 

We recognise the fundamental BCJ relation for (n — k + 2) external particles. For k > 3 we have 
(n — k + 2) < n. We may therefore use the induction hypothesis and we conclude 

Ri = 0. (83) 


The argument for R n is very similar. We first exchange the summation indices as 

n—In —2 n—2 k 

EE/(^) = EE/(^)’ 

i=2 /:=/ k=2 i=2 


(84) 


We then obtain 


n—2 k / i 

Rn = I ^P2Pj )A n (i,3,...,i,2 g ,i+l,...,k,P\-P,k+l,...,n-l,n) 


k=2 7=2 \;=1 


= 0. 


(85) 


Here we used momentum conservation in the sum over j. Again we recognise the fundamental 
BCJ relation in the form of eq. ©. It follows that R n vanishes. 

Exchanging the summation indices for R 2 one obtains 


n—2 n —1 

* = E E 

k=2 I=k+1 


( 86 ) 


l / n 


X L L 2 P 2 Pj )A n (k+l,...,i,l g ,i+l,...,l,P\-P,l + l,...,n-l,n,l,3,...,k). 

i=k \j=i +1 


We may split the sum over j as 


\ 1 r\ a. a \ 1 r\ a a . \ 1 /\ /v 

L 2 P2Pj = l 2 PiPj + L 2 P 2 Pj 

j=i+\ j=i +1 J=/+1 

A B 

The terms of type A vanish again by the induction hypothesis 

Ilf L 2 P 2 PJ ) A i-k+2(P,k+l,...,i,5.g,i+l,...,l) = 0. 
i=k \j=i +1 / 


(87) 


( 88 ) 


Note that the sum over i extends only to (/ — 1), the case i — l contributes only to the terms of 
type B. 
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For the terms of type B the sum over j is independent of i and may be taken outside the sum 
over i. The sum over i vanishes then due to the U (1 j-decoupling relation, given in eq. ©: 


( £ 2 P*Pj J L A ^-*+2(A^+l,-,b2 g ,/+l,= 0. (89) 

\j=l +1 / i=k 

We therefore conclude that 


R 2 = 0. (90) 

Putting the partial results for R\, R 2 and R„ together we find that 

4(0) = 0. (91) 


This completes the proof. 


6 Conclusions 

In this paper we provided a proof of the fundamental BCJ-relation, stated in eq. ©, for primitive 
tree amplitudes in QCD. The proof holds for massless and massive quarks. For the proof we used 
induction and BCFW-recursion relations. 


A Reference spinors for massive particles 

A.l The helicity configuration l + ,2+,rc~ 

In this appendix we show that the system 

KilgiT) = yi p b n+y K„(g„+| = yi (p\+ +y 2 {p 2 +\ 
has a solution. Expressing (p 2 + \ in terms of (l\ + \ and ( l n + \ 

(P2 +1 = £#('■+l + M<t.+l, 


[44 


[44 


we obtain the system of equations 


Ki = y\c n , 

KiXt — y\C n QL\h ni 

[p 2 ln 


= y\c\ +y 2 - 


[44] ’ 


^nKn yicicx /z A.i Ty2 


[hp 2 \ 
[44]' 


(92) 


(93) 


(94) 
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The variables ai and 0 C 2 are defined in eq. (l26l) . the variables c 1 and c„ are defined in eq. (l32l) . 
We look for a solution for the variables Ki, k„, y\, y 2 , An and A„. A possible solution is 


Kl 

5 

C 1 

1 

K„ 

yt 

1 

C 1 

J2 


- 1 + 


2 P2 ln 


2hl n ’ 

( Plln) 

(hln) ’ 


(95) 


Al = 


Pn (h + I// 2 I ln+) 


A„ 


2l\l n ( l\ + \p/l \ /«+) 


(2hl n ) 2 -p 2 lP l + (2hl n )(2p 2 l n y 

A.2 The helicity configuration 1+2" n 

In this appendix we show that the system 

Kiki+) = yi\P 2 +)+y n P b n + 


{2l\l n ) — PiP„ + (2/i/„) (2p 2 l n ) 


K n (q n +\ — y n (p\+ 


has a solution. Expressing \p 2 +) in terms of |/i+) and |/„+) 

Ib+> = 


</l/n) 


(W 


we obtain the system of equations 


K„ = 

K77A77 


y n c \ QL n '\j \, 

(hpi) 


Kl = ynCn+yi 


K 1 X 1 = 


(hln) ’ 


A possible solution is 


Ki = 1 + 


Ai 


J2 = 

2 / 1/77 (/l + |i/2| /;;+) 


y n c n U\ 7-tj T" y2 

(i>2/«) 

( /1 //7 ) 

2/l i?2 

Cl 

2/,/„’ K " “ 

? 

Cn 

[/lP2] 

1 

[/1 //7 ] ’ 

5 

C 77 


(96) 


(97) 


(98) 


(99) 


(2hl n y-p\p 2 + {2hln)(2hp 2 ) 


A ,7 — 


P\ (h + |i/2| //?+) 


(2l\l n ) — P\Pn + (2/i/n) (2l\p 2 ) 


B The large z-behaviour in the eikonal approximation 

Let us consider a theory with massless or massive scalar “hard” particles, denoted by a hat and 
QCD-like “soft” particles (gluons, quarks, anti-quarks), denoted without a hat. The momenta of 
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the hard particles are of order z 1 , the momenta of the soft particles are of order z°- The Feynman 
rules for this toy theory are as follows: The hard particles interact only through three-valent 
vertices. The Feynman rule for the three-valent vertex involving three hard particles with the 

/v/v/v 

cyclic order (1,2,3) is simply i, for the cyclic order (1,3,2) we have (— i). Furthermore there is 
a three-valent vertex, involving two hard particles and one soft gluon. The Feynman rule for the 
cyclic order (1,2,3) reads 


1 (/ 3 i - P3) ■ 


( 100 ) 


There are no vertices involving only one hard particle. The Feynman rules for the vertices in¬ 
volving only soft particles are the standard (cyclic-ordered) QCD Feynman rules, listed in ap¬ 
pendix o 

Let us consider the situation of three hard particles 1, 2 and h and (n — 3) soft particles 3, 
..., (n — 1). We assume particle 2 to be massless and particles 1 and h to have the same mass m 
(which may be zero). We will denote an amplitude in this toy theory by 


1 eikonal 


(i,3,...,j,2,/+l,...,n—l,n), 


( 101 ) 


and we define 


reikonal 


(z) 


£ \2p2Pn+ £ lP2Pj W konal (t,3,. 
i =2 \ j=i +1 / 


1 , 2 , i+ 1 , ...,n — l,n). ( 102 ) 


We would like to show that 7® lkonal (z) goes to a constant for large z. Then the quantity 

O 3 7,f konal (z) (103) 

with O 3 defined as in eq. (l69l) falls off like 1 fi¬ 
ll will be convenient to introduce soft currents 

•Ct(a,-,*0, (104) 

involving {b — a+ 1) soft on-shell particles a,a+\,...,b and one soft off-shell gluon leg. The 
momentum of this soft current is 


P 


E Pk- 


k=a 


(105) 


We may group the Feynman diagrams contributing to 7® lkonal (z) into sets, where exactly r soft 
currents couple to the hard particles 1, 2 and h with 1 < r < n — 3. Therefore we have a decom¬ 
position 


7,f onal ( Z ) = £ 7 eik° nal ( z )_ (106) 

r= 1 
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2 

«iQQz Jf l 
J s a f _i 


n 


Figure 1: Diagrams showing how the two outermost soft currents 7^ oft and ma y couple to 
the hard particle 2. There may be further soft currents attached to the hard lines below the dashed 
line. 


We will show that each contribution 7® lkonal (z) individually goes to a constant for large z. 

Let us discuss 7 a ', konal (z) with r soft currents 7j oft , 7* oft and associated momenta P \, P r . 
The cyclic order among the soft currents is respected in each diagram contributing to 7“ konal (z). 
We will use the notation 


b 

Pa,a+i,...,b Pk- 

k=a 


(107) 


Let us first discuss the situation, where two or more soft currents couple to the hard line 2. These 
contributions add up to zero in 7® lkonal (z). In order to see this, consider the situation, where the 
two outermost soft currents coupling to 2 are 7* oft and There are four possibilites, how 
these soft currents may couple to 2, shown in fig. ([[]). With the inclusion of the prefactors of the 
fundamental BCJ-relation, these contributions add up to zero. 

Let us now consider the case, where one soft current 7* oft couples to the hard particle 2. Here 
we get the contribution 


2 do / sokt 

~(2p 2 P a ) P ~ a = 2p l J s / t + 2p n r/ i +0(z°). (108) 

(P2+Pa) 

We may now add up all contributions and obtain 


reikonal 

1 n,r 


(z) 


L ( 2 p 2 P” + E 

i =0 \ j=i +1 




i= 1 




(-i)' 



2pVf ft \ 

2piPi,...,k J 



wrM 
2piPi,...,k J 

2/3,;7 / solt \ 

2 PnPl,...,r J 


+0(z°). 



2p n Pl,...,r J 


(2p 1 7 ; soft + 2p ;; 7f oft ) 


(109) 


The terms in the first sum come from diagrams, where all soft currents couple either to the hard 
particle 1 or n, the terms of the second sum correspond to diagrams, where exactly one soft 
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current couples to the hard particle 2. Noting that 


2 P2Pn+ £ 2 p 2 Pj = 2p]P] j — 2p n P i+ ] j-+ O (z°) (110) 

j=i +1 


one sees that 


reikonal 


(z) 



( 111 ) 


as claimed. 


C Cyclic-ordered Feynman rules 

In this appendix we give a list of the cyclic-ordered Feynman rules. They are obtained from the 
standard Feynman rules by extracting from each formula the coupling constant and the colour 
part. The propagators for quark and gluon particles are given by 


= 

The cyclic-ordered Feynman rules for the three-gluon and the four-gluon vertices are 


. pj + m 

l T T 1 
p - — m z 

-ig^ 


( 112 ) 



i (pA'B -p^) + g W 3 +g fJ w -p^)] : 



(113) 


The Feynman rule for the quark-gluon vertex is given by 



( 114 ) 
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